Probe and Control of the Reservoir Density of States in Single-Electron Devices 
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We present a systematic study of quasi-one-dimensional density of states (DOS) in electron ac- 
cumulation layers near a Si-Si02 interface. In the experiments we have employed two conceptually 
different objects to probe DOS, namely, a phosphorus donor and a quantum dot, both operating 
in the single-electron tunneling regime. We demonstrate how the peaks in DOS can be moved in 
the transport window independently of the other device properties, and in agreement with the the- 
oretical analysis. This method introduces a fast and convenient way of identifying excited states in 
these emerging nanostructures. 

PACS numbers: 73.21.-b,61.72.uj, 83.35.-p 



The appearance of discrete energy spectra is often re- 
garded as a fundamental property of quantum systems. 
However, the limit of vanishing energy spacing can be 
met in many mesoscopic systems at the interface of the 
quantum and classical regimes meaning that the energy 
levels can be treated as a continuum. This gives rise to 
the concept of density of states (DOS), the integral of 
which over an energy domain yields the number of quan- 
tum states in that region. The DOS concept has been 
successfully applied in explaining transport, absorption 
and emission, and quantum statistical phenomena in such 
devices [H. In particular, DOS is a key element in the 
low-temperature behavior of metal-oxide-semiconductor 
field-effect transistors [H (MOSFETs) which constitute 
the cornerstone of information processing circuits today. 
Since device miniaturization has now reached the point 
where quantum effects of single atoms can dominate the 
operation characteristics [3, 0, S S| , there is an urgent 
need to understand and distinguish continuum DOS ef- 
fects from the discrete quantum behavior. In this paper 
we present a systematic study of the reservoir DOS in 
two gated MOSFET nanostructures with quantum chan- 
nels defined by either a quantum dot or the extreme case 
of a single donor atom. 

In a silicon MOSFET a positive gate voltage is ap- 
plied to induce an electron layer directly below the Si- 
Si02 interface. Hence, the electron dynamics is essen- 
tially limited to two dimensions, leading ideally to a con- 
stant DOS. In a narrow channel, however, the contin- 
uum approximation is only valid in the longest direc- 
tion, giving rise to quasi-one-dimensional (QID) DOS 
with highly non- uniform characteristics, see Fig. [TJa). 
Conductance modulations attributed to QID density of 
states were first observed in a parallel configuration of 250 



narrow MOSFET channels and later in a single chan- 
nel d, [^, . Since the cumulative conductance through 
all occupied states was measured in these experiments, 
the studies were limited to rather low electron densities 
in the channel. The QID density of states is also in- 
timately related to ballistic electron transport through 
quantum point contacts 11, IT^ . and has been probed 



in single carbon nanotubes using scanning tunneling mi- 
croscopy 
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FIG. 1: (Color online) (a) Schematics of single-electron tun- 
neling through a discrete quantum level, /x_d, together with 
the reservoir Fermi levels, Hl/r, and density of states, (b) 
Schematic stability diagram showing conductance through the 
discrete quantum level as a function of the source-drain bias 
and the plunger gate voltage. The numbered lines of increased 
conductance arise from the corresponding DOS peaks in panel 
(a), (c) Zeeman splitting of the reservoir and discrete energy 
levels in a magnetic field. 



In single-electron transport through discrete quantum 
states m , the current is directly proportional to the reser- 
voir DOS at a given energy, see Fig. [TJ Considerable ef- 
fort has been directed towards the study of local density 
of states of a reservoir in the vicinit y of an impurity atom 



in a GaAs quantum well [15|, [16|, lUl- Here, the local 



DOS was dominated by disorder due to impurity scat- 
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tering, resulting in reproducible but irregular features in 
the conductance which behaved in a complicated way as 
a function of magnetic field. Thus the term local DOS 
fluctuations was introduced. In our case, the behavior is 
less complicated, and hence we do not refer to peaks in 
DOS as fluctuations. 

In gated devices, local DOS effects have also been ob- 
served in tunneling through a vertical quantum dot [l^ , 
for which a schematic stability diagram is shown in 
Fig. [Ijb). Recently, gated quantum dots have attracted 
great interest due to their tunability, and lines in the 
stability diagrams not attributed to excited states have 



been observed in various structures [19|, |20|, l2l|. In 
Refs. [13, 21 1 however, these lines were due to phonon 
modes, whereas in Ref. [l3| they were suggested to arise 
from reservoir DOS but were not studied thoroughly. 

Motivated by the Kane proposal [23] for a quantum 
computer based on shallow donors in silicon [23|, there 
has been considerable development in electron transport 
through gated single donors 0, 0, Q ■ In these exper- 
iments, conductance lines attributed to reservoir DOS 
have been observed but, again, not studied in detail. In- 
deed, it has not yet been proven that these features are 
due to the reservoir DOS. In this paper we present a 
systematic study of the reservoir DOS in two different 
nanostructures, namely, in a recently introduced double- 
gated single-donor transistor (Structure A) and in a 
novel multi-gated silicon quantum dot ^231 (Structure B), 
shown in Fig. [2] In contrast to previous studies of DOS, 
the double- and multi-gated designs allow us to map con- 
veniently the reservoir DOS at a wide range of energies 
and electron densities. 
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FIG. 2: (Color online) (a) Top view of Structure A and (b) 
Structure B and a schematic cross section of (c) Structure A 
and (d) Structure B in the y = plane. The insets show 
schematic potential landscapes along x. 

The high-purity [001] silicon substrates for both device 
structures went through similar fabrication steps to those 
described in detail in Ref. 0. In Structure B, however, 
there are three layers of metallic gates deposited with no 
implanted donors underneath. The working principle of 



both devices is that an accumulation layer of electrons 
is induced directly below the Si-Si02 interface using the 
lead gate voltage V\g. This layer constitutes the reservoirs 
for single-electron tunneling through a single donor or a 
quantum dot, the electrochemical potentials of which can 
be tuned by the plunger gate voltage V^g. For Structure 

A, the plunger gate works also as a barrier gate, depleting 
the reservoirs in the vicinity of the donor. For Structure 

B, however, we have separate barrier gates which can be 
used to tune the coupling of the quantum dot to the left 
and right reservoir independently. The electron densities 
in the reservoirs can also be controlled independently, but 
for simplicity, we keep them at the same value here. As 
the plunger gate voltage is increased, the electrochemical 
potential of the donor or dot, iiu, shifts down, eventually 
entering the source-drain bias window leading to single- 
electron tunnehng through the device, see Fig. [TJa-b). 

The sequential tunneling rate is directly proportional 
to the reservoir density of states, which leads to a peak 
in the source-drain current if fio is aligned with a peak 
in the density of states. Critically in our devices, we can 
shift the Fermi levels of the reservoirs with respect to 
the conduction band minima by changing the lead gate 
voltage, T^g, which in turn moves the DOS peaks with 
respect to the transport window. The effect of the lead 
gate voltage on /x^ can be compensated by the plunger 
gate [2^. Therefore, if we observe a conductance peak 
to shift with Vig in the stability diagram, we can iden- 
tify it to arise from a peak in the reservoir DOS, clearly 
distinct from features due to excited states in the dot 
or donor which do not move with Vig. Previously, DOS 
peaks have been probed by changing the temperature |l3] 
or magnetic field [2^ , both of which also have an impact 
on the donor or the dot, and which are many orders of 
magnitude slower than the method presented here. 

The measurements were carried out in '^He-'^He dilu- 
tion refrigerators below 100 mK temperatures. The con- 
ductance was measured using standard lock-in techniques 
and the direct current was measured from the same sig- 
nal after low-pass filtering the modulation arising from 
the 10-50 fiV lock-in excitation. For Structure B, only 
the direct current was recorded and the differential con- 
ductance was extracted numerically. 

Figure [3] shows measured stability diagrams for Struc- 
ture A [panel (a)] and B [panel (d)]. In Structure A 
the coupling of the donor to the left reservoir is much 
weaker than to the right reservoir, and hence only lines 
with positive slopes corresponding to the left reservoir 
are visible. This is justified by the sequential tunneling 
model, in which the current through the device is given by 
/ = 2erinrout/(2rin + rout), where Fin/out is the in/out- 
tunneling rate. Thus the current is determined by the 
small rate, independent of whether it corresponds to in- 
or out-tunneling. The tunneling rates in Structure B can 
be tuned to be symmetric, but considerable asymmetry 
persists in the working point of Fig. [3] 
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response of the Fermi level Ep to the lead gate voltage 
to be A{Ef - £^c)/AMg = (11-3 ± 1.5) meV/V, where 
Ec is the conduction band minimum. We have modeled 
this response using technology computer aided design (27l | 
taking into account the three-dimensional geometry of 
the system. By matching the measured response with 
the model, we obtain an estimate of the absolute value 
of the Fermi level to be Ep — Ec = (25 ± 5) meV. 

The average energy spacing between the DOS peaks is 
obtained from Fig. [3] to be (0.33 ± 0.03) meV for Struc- 
tures A and B. We model this by calculating the single- 
particle energy spectrum for a triangular well Vz{z) = 
eF\z\/Q{—z) in the z direction {Q is the Heaviside step 
function), an infinite W = 70 nm box potential in the y 
direction, and a semi-infinite box potential in the x direc- 
tion. Here, F is the electric field in the vicinity of the Si- 
Si02 interface due to the lead gate. The one-dimensional 
DOS in the x direction, for which the continuum approx- 
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imation is valid, is given by Di^{E) = j {j^hy E\ 

where m^, is the effective mass of the electron in the x 
direction. Thus we obtain for the total density of states 

E>\.E) = — > ^ , , , (1) 

^ ^E - Ey,z{ny,nz,n^) 



riy ,nz .n^ 



FIG. 3: (Color online) Conductance of (a) Structure A with 
Vig = 1.87 V and By = 4 T, and (d) Structure B with VIg = 
3.00 V, Vib = 0.52 V, Vrh = 0.68 V, and By = T as a function 
of the source-drain bias and the plunger gate voltage. The 
dashed line in panel (a) represents the trace shown in panels 
(b) and (c) for l^d = —2.5 mV as a function of the lead 
gate voltage for differential conductance and direct current, 
respectively. Panels (e) and (f ) correspond to the dashed trace 
shown in panel (d) for AVpg = 2.1 mV from the transport 
window at Vsd ~ 0. The capacitive coupling of the lead gate 
to the donor or the dot has been compensated by the plunger 
gate in the traces. Panel (g) shows the theoretical prediction 
for DOS corresponding to panels (b,c). Here, Eq. has 
been convolved with the Gaussian exp{—Ep /[2{500 /ieV)^]} 
to account for broadening due to experimental nonidealities. 



The stability diagrams shown in Fig.[3]Ja,d) do not pro- 
vide enough information to distinguish whether a spe- 
cific conductance feature is due to the reservoir DOS 
or an internal excited state. However, the lead-gate- 
compensated traces shown in Fig. [3l^b,c,e,f) reveal the 
origin of the conductance peaks. Most of the peaks 
move with the lead gate voltage, and hence correspond to 
the reservoir DOS. However, Fig.[3fb,c) shows a feature 
0.4 meV above the ground state for all Vig. As justified 
by measurements carried out in various magnetic fields 
(see Fig. HI, this line corresponds to the spin excited state 
of the neutral donor orbital ground state. Similarly, the 
vertical lines in Fig.[3l[e,f) are attributed to excited states 
of the dot. 

For Structure A, the slope of the DOS lines yields the 



where rij,, n^, and Uy are the indices for the energy levels 
in the y, z, and valley degrees of freedom, respectively. 
The six-fold valley degeneracy of bulk silicon is lifted due 
to different effective masses nit = 0.190 x toq and m/ = 
0.916 X Too for the transverse and longitudinal directions, 
respectively Here, TOq is the electron rest mass. Thus 
the effective mass in the x direction depends on the valley 
degree of freedom. According to Eq. ([1]), a peak in DOS is 
observed at energies matching the single-particle energy 
levels in the two-dimensional yz potential, i.e., for E — 

The peak spacing resulting from the confinement in 
the y direction is given by Tr^h?{ny + l/2)/{myW^). 
This yields the minumum peak spacing of 1.2 meV for 
the transverse and 0.25 meV for the longitudal effective 
masses. Each excited state in the valley and z degree of 
freedom superimposes an additional series of peaks, and 
hence the resulting DOS structure can be rather compli- 
cated. Figure [3]^g) shows the DOS in the vicinity of the 
Fermi level obtained above for Structure A. The electric 
field of is chosen to be F = 2 MV/m to match the peak 
separation observed in the experiments. This field is ap- 
proximately an order of magnitude lower than expected 
due to the crude model used to calculate the peak spac- 
ing. This simple model does not include the true details 
of the trapping potential including disorder or the elec- 
tric field induced by the accumulation layer, which calls 
for deeper theoretical analysis. Furthermore, we have not 
included the possible splitting of all valley degeneracies 
which would further reduce the average peak spacing and 
consequently, increase the required electric field. 
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Figure[TJc) shows a schematic diagram of how the DOS 
peaks behave in magnetic field B = By. The Zeeman ef- 
fect shifts the kinetic energy at the Fermi level {Ep — Ec) 
for the spin-down electrons up by and vice 

versa for spin-up electrons. Here, is the electron g 
factor in the reservoirs and iib = eh/{2mo). Thus the 
effect of this term is to shift and split the DOS peaks. 
However, if the electron g factor for the donor or dot, 
go, is almost the same as gL, this compensates for the 
splitting and results only in a downwards shift of the 
DOS peaks by gL^BB/2 with respect to the transport 
window. In our case, the DOS peaks simply shift and 
do not split in magnetic field, as demonstrated in Fig. HI 
which implies an ~ gL ~ 2 consistent with previous mea- 
surements [28|. In contrast, Zeeman splitting of DOS 
peaks has been observed in tunneling through a Si shal- 
low donor in a GaAs/AlAs/GaAs junction [26i |. 

Figure lUJd) shows the energy shift 5E of the conduc- 
tance peaks with respect to the transport window as 
a function of magnetic field. We observe that within 
the estimated error, the shift of the DOS peaks and the 
spin excited state is in agreement with 5E = fisB and 
5E = 2^bB, respectively. Thus our observations can be 
explained by the Zeeman effect alone. 
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FIG. 4: (Color online) (a-c) Stability diagrams of Structure 
A in magnetic fields ranging from 1 to 6 T . The spin excited 
state [L6 in panel (a)] is marked with a white and blue arrow, 
and the DOS lines are labeled with green arrows, (d) Shift 
of the conductance lines in magnetic field. Data for each line 
is ofi'set in energy for clarity. The dashed lines show a slope 
2^B for L6 and /xb for the other data. 

In conclusion, we have reported the first studies of 
reservoir DOS in nanostructures through gate control 
of the reservoir electron density. Magnetic field spec- 
troscopy of the DOS peaks revealed behavior consistent 



with energy shifts due to the Zeeman effect alone. Both 
the single donor and the quantum dot devices studied 
exhibited similar conductance patterns arising from the 
reservoir DOS. Our findings not only confirm the inter- 
pretation that additional conductance lines observed pre- 
viously in transport through single donors are due to 
reservoir states [1, 0, Q j but also demonstrate that these 
states can be controlled and studied in a quantitative 
manner. 
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